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Abstract 

We study the critical properties of the weakly disordered two-dimensional Ising 
and Baxter models in terms of the renormalization group (RG) theory general- 
ized to take into account the replica symmetry breaking (RSB) effects. Recently it 
has been shown that the traditional replica-symmetric RG flows in the dimensions 
D = 4— e are unstable with respect to the RSB potentials and a new spin-glass type 
critical phenomena has been discovered |11],|12]. In contrast, here it is demon- 
strated that in the considered two-dimensional systems the renormalization-group 
flows are stable with respect to the RSB modes. It is shown that the solution of 
the renormalization group equations with arbitrary starting RSB coupling matrix 
exhibits asymptotic approach to the traditional replica-symmetric ones. Thus, in 
the leading order the non-perturbative RSB degrees of freedom does not effect the 
critical phenomena in the two-dimensional weakly disordered Ising and Baxter 
models studied earlier. 
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The effects produced by wealc quenched disorder on the critical phenomena near the 
phase transition point have been studied since many years ago According to the 

Harris criterion the disorder effects the critical behaviour only if a, the specific heat 
exponent of the pure system, is positive. In this case a new universal critical behavior, 
with new critical exponents, is established sufficiently close to the phase transition point 
for (T/Tc — 1) = r <^ r„ = u^/" 0, where u -C 1 is the parameter which describes the 
strength of the disorder. In contrast, when a < 0, the disorder appears to be irrelevant 
for the critical behavior. 

Originally the modified critical behaviour has been derived for the classical 0^ model 
near four dimensions , and later it has been studied for the two-dimensional Ising , 
Baxter [Q and Potts |^ models by various renormalization group (RG) techniques, and 
by numerical simulations 0. 

In dealing with the quenched disorder the traditional approach is the replica method, 
and in terms of replicas all the results obtained for the systems listed above correspond 
to the so-called replica-symmetric (RS) solutions. Physically it means that olny unique 
ground state is assumed to be relevant for the observable thermodynamics. The prob- 
lem, however, is that in the presence of the quenched disorder there exist numerous 
local minimum configurations separated by finite barriers, and in this case the direct 
application of the traditional replica-symmetric RG scheme may be questioned. 

On the other hand, it is the Parisi Replica Symmetry Breaking (RSB) scheme which 
has been developed specifically for dealing with disordered systems which exhibit nu- 
merous local minima states (see e.g. [0). Recent studies show that besides mean-field 
theory of spin-glasses the RSB approach can also be generalized for situations where one 
has to deal with fluctuations as well p[,0,|jlO[- In the paper |]ll| qualitative arguments 



were presented demonstrating the mechanism how the summation over multiple local 
minima configurations could provide additional non-trivial RSB interaction potentials 
for the fluctuating fields. 

One can hope that such type of the generalized RG scheme self-consistently takes 
into account relevant degrees of freedom coming from the numerous local minima, and in 
particular, the instability of the RS fixed point with respect to the RSB would indicate 
that the multiplicity of the local minima is relevant for the critical properties in the 
fluctuation region. And inversely, if the traditional RS RG flows turn out to be stable 
with respect to the RSB modes, then it can be concluded that such non-perturbative 
degrees of freedom are irrelevant for the critical behaviour in a considered system. 

The first example of replica symmetry broken solutions in the renormalization group 



has been suggested in |]13[ in the context of the 2D random field XY-model, where 
instability of the RG flows with respect to the RSB modes has been discovered. Similar 
phenomena has been observed in the weakly disordered classical ferromagnet (the 0^- 
theory) in the dimension D = 4 — e Here the RSB degrees of freedom produce 

dramatic effect on the asymptotic behaviour of the RG flows, such that for a general 
type of the RSB there exist no stable fixed points, and the RG arrives into the strong 
coupling regime at the finite scale (~ exp(— [jl2|. 



On the other hand, similar considerations for weakly disordered 2D Potts model 
shows that, although the traditional RS fixed point also turns out to be unstable, there 



2 



exist stable non-trivial fixed point characterized by the continuous RSB structure of the 



coupling matrix |14 



In this Letter we report on the RSB solution for the two-dimensional Ising and 
Baxter models with random bonds. We consider the models with spin couplings having 
small fluctuations around a mean ferromagnetic value. This gives a possibility to study 
the model in the continuum limit, because one approaches the critical point sufficiently 
close before the randomness becomes relevant. For the two-dimensional Ising and Baxter 
models this allows to use the renormalization group based on the fermion representation. 

In contrast to the previous studies, in the considered systems the KG flows turns out 
to be stable with respect to RSB modes. The explicit solution of the corresponding RG 
equations shows that for any starting RSB structure of the coupling matrix the asymp- 
totic solutions become replica-symmetric. It means that the RSB (" non-perturbative" ) 
degrees of freedom appear to be irrelevant for the critical phenomena in these systems, 
and they exhibit the usual replica symmetric critical behavior which has been studied 
earlier ||,||. 

The Ising Model 

It is well known that the two-dimensional Ising model in the critical region is equiv- 
alent to the problem of the free fermions [|T^. The effect of the quenched disorder can 



be described by the random contribution to the effective temperature r. Thus, the ran- 
dom bond Ising model in the critical region can be reduced to the field theory with the 
following Lagrangian: 



d^x[^p^^p + {t + 5t)^/j^], (1) 



where 6t is the Gaussian random variable with (5t^) = 2u <^ 1, and ^p and ^p = iipay 
are the two-component real fermion fields. 

The averaging over quenched random variable 5t can be performed in terms of the 
standard replica procedure. The resulting replica Lagrangian has the following form: 

= -j £x[- J2rid+ T)r - T E ga^rr^V], (2) 

a=l a,b=l 

In the usual replica-symmetric theory the coupling matrix does not depend on the replica 
indices: Qab = u. Then, in terms of the standard RG scheme one can easily calculate the 
critical properties of the system, and in particular for leading singularity of the specific 
heat one finds: C(r) ~ loglog(l/|r|) §. 

According to the discussion in the introductory part of this Letter, in our present 
scenario we assume that due to non-perturbative effects the replica symmetry in the 
coupling matrix Qab can be broken. In other words, the starting point of the further 
analysis is the assumption that the matrix gab has the Parisi-type structure (see e.g. 

i)- 

The corresponding (one-loop) RG equations for the replica matrix g^b and for the 
mass parameter r are: 
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= - 9ac9cb- (3) 

dlogT 1 ^ 

where ^ is the standard RG rescahng parameter. 

In terms of the Parisi RSB scheme in the hmit n the matrix Qa^b is 
parametrized by the function g{x) defined in the interval < x < 1 (replica-symmetric 
situation corresponds to g{x) = const, independent of x). According to the standard 
technique of the Parisi RSB algebra (see also |[Tl|] , [|12| ) the equation can be rep- 
resented as follows: 

"^^^""^ - \[xg'ix) + 2g{x) f dyg{y) + f dyg\y)], (5) 



d£^ TT Jx Jo 

The structure of this equation is similar to that for the 0^-theory with the broken 
replica symmetry, and it can be solved following the method suggested in [jl2|. Simple 
calculations yield: 

>( .^ ^ g'{x,0)exp{-^J^d7]g{7])) 

g{0,O= 9(0,0) exp{-- g{v)dv)- (7) 

TT Jo 

where g'{x,^) = £g{x,0 and g{^) = Jq g{x,^)dx. 

Equations (^-(0) allow to study the large scale asymptotic behaviour of the renor- 
malized function g{x,^) provided the starting function is given. Simple calculations 
show that at any value of x in the limit ^ ^ oo the derivative in (P) tends to zero 
as ['C(logO^]^^! while g{^) ~ It means that the asymptotic renormalized function 
g{x,^ oo) tends to become flat in the whole interval < x < 1). In other words, for 
any starting RSB function g{x, ^ = 0) the asymptotic solution of the RG eq.(|^) becomes 
replica-symmetric. 

One can also easily check that slowly decaying corrections to the RS solution does not 
effect the leading singularity of the specific heat. From the eq.(P for the renormalized 
temperature parameter one finds: 

r(O = r(0)exp(-- g{v)dv)- (8) 

TT Jo 

Then, using eqs. (P),(0) and (§) one can derive the following closed form relation on 



r(0) ^(0,0) f^(v).2^n( a J^^) (a\ 
In the asymptotic limit ,^ ^ oo one finally gets: 
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r(0~4 + O(^=^). (10) 

which in the leading order coincides with the RS resuh. Correspondingly, the leading 
singularity of the specific heat (which can be estimated as C(r) ~ f""^^^^^^^ d^r^^) 
remains the same as in the RS case as well. 

The Baxter Model 



The Baxter model [16] can be formulated in terms of two 2D Ising models coupled 
by four-spin interactions. The strength of this coupling is described by parameter A 
(the case of A = corresponds to two independent 2D Ising models). The pure Baxter 
model can be solved exactly and in particular, its specific heat exponent is proportional 

to A: CpureiT) ~ |t| (for |A| <^ 1). 

Similar to the weakly disordered 2D Ising system considered above, in the continuum 
limit near the critical point the corresponding Baxter model can be described in terms 
of the two-component complex fermion field: 

m = -J d'xi^^pd^ + i(r + 5r(x))# - iA(#)(^^)] (11) 

After averaging over the random function 6t{x) one gets the following replica La- 
grangian: 

a=l a,b 

In the RS case for the replica coupling matrix one has: gab = A^ab + u- The critical 
properties of such system has been studied earlier , and it was shown that for A > 
(when the specific heat of the pure system is divergent) the leading singularity of the 
specific heat coincides with that of the random Ising model: C(r) ~ Inln^. On the 
other hand, for A < (when the specific heat of the pure system is finite), the singularity 
of the specific heat changes to: C(r) ~ |r| ^, where A* = A exp(— m/| A|). 

Here we study the situation when the coupling matrix Qat has a general RSB Parisi 
structure. The corresponding RG equations are: 

= --{gabigaa + 9bb) - 2j29ac9cb} (13) 
d^ TT V 

dilnr) 1 

According to the standard technique of the Parisi RSB algebra, in the limit n ^ 
the matrix Qab is parametrized in terms of its diagonal element g and the off-diagonal 
function g{x). The replica symmetric situation corresponds to the case g{x) = const 
independent of x. After simple algebra, instead of the eqs.(|T3| ) , (|T^) in the limit n — > 
one gets: 
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P = -I lo' dxg^x) 

igix) = Bg{x) - 2gg{x) - dy[g{y) - g{x)]'] (15) 

|(lnr) = 1(^-2^) 
where g = Jq dxg{x). 

The solution of the eqs. ([ISD can be represented in the same way as in the Ising case 
considered above (cf eq.(p|)): 



^ g'jx, 0) exp{-| dri{2g{r]) - g{r]))} 

^^"^ ' [i + IIoy9'iy,o)dyj,^dvexp{-lJSde{2m-m)}?' ' 

Using the above equation one can easily study the asymptotic (for ^ oo) behaviour 
of the function g{x,C,)- The behaviour of the solution depends on the sign of the starting 
parameter g^ = g{^ = 0) — (^(^ = 0) which corresponds to the initial coupling A in the 
replica symmetric case. Consider the two cases separately. 

1) ^70 > 

The corresponding replica-symmetric RG equations have the following asymptotic solu- 
tion ^■. 

2 r 1 111 a 



QRS — ttL^ glngJ ' 5(ln5)2 

(17) 

s 2 rl 2_i I b 

gns — gingJ 5(111^)2 
where a and b are constants. Then, using eq.(|l6D for the RSB correction to these 
solutions: g{x,^) = gnsiO + 9RSBix,0 and g{^) = gnsiO one immediately 

finds that 

(18) 

where V{x) is a function which is defined by the starting RSB function g{x, ^ = 0). This 
means that on large scales the RSB deviations are vanishing in comparison with the RS 
solution, and the asymptotical properties of the model are not affected by the starting 
RSB. 

2)go<0 

In this case the replica-symmetric RG trajectories arrive to the fixed point {g^,, g^,), which 
depends on the initial coupling parameters g{C, = 0) = g{0) and go [Q: 



g* = -|5'o|exp 



1 90 



(19) 



Using eqs. ([T5|) or (|l^) this fixed point can be easily shown to be also stable with respect 
to general RSB deviations. In other words, the RG trajectories defined by eqs.(|15|) arrive 
to the same fixed point {g^, 0) as in the RS case. 
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In conclusion, we have studied the effects of the rephca symmetry breaking in the 
renormahzation group for the 2D random Ising and Baxter models. We have found that 
the traditional replica-symmetric renormahzation group flows are stable with respect to 



the RSB modes. Thus, unlike the other suggestions |Tj, we conclude that the effects of 
possible non-perturbative degrees of freedom in the 2D random Ising and Baxter models 
must be irrelevant for the critical behaviour studied earlier. 
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